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Abstract. - Corrections to tlie atomic orbital sum rule for circular magnetic x-ray dichroism 
in solids are derived using orthonormal LMTOs as a single-particle basis for electron band 
states. 



Atomic physics affords a theory of x-ray dichroism by providing a set of sum rules which 
relate dichroic intensities, integrated over a finite energy interval, to the ground-state expec- 
tation value of effective one-electron operators For circular magnetic x-ray dichroism, 
that is the difference in absorption between right- and Icft-circularly polarised photons in 
a system with a net magnetisation, the effective operators coincide with spin and orbital 
multipoles 

The spherical symmetry and the discreteness of the spectrum governing the atomic results 
do not hold for an atom in a solid, where the spin and angular-momentum character selected 
by a particular x-ray transition is spread out over a band of final states. This difference 
has hindered the identification of a well-defined connection between the atomic sum rules 
and band-structure calculations of magnetic x-ray dichroism Q. Except in cases of strong 
electronic correlations, such calculations have been very successful in simulating experimental 
absorption spectra (^,0. It therefore seems desirable to derive a band-structure formalism 
which exhibits the atomic sum rules as the dominant term. This should be important not 
only for the interpretation of x-ray dichroism in cases where current density-functional band 
theory works, but also as a prerequisite for understanding dichroism in strongly correlated 
materials. The current paper is an attempt in this direction. 

By leaving a localised hole, inner-shell photo absorption selects a specific site in the solid, 
which we shall label by i? = 0. A local process is thus expected to control the excitation, 
to leading order. Additional contributions should emerge when the remaining sites in the 
lattice are taken into account, that is, when electron delocalisation is included. In this case, 
a minimal set of orthonormal linear muffin-tin orbitals (LMTOs) provides a suitable single- 
particle basis P-p^. 
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The macroscopic quantity of interest is the polarisation and energy dependent extinction 
coefficient, k^{uj) JT^ . In a microscopic description, this is given by 

where f^{u>) stands for the forward scattering amplitude, as determined by the p • A coupling 
between x-rays and electrons. Photon energy (in units of Ti) and polarisation are identified 
by to and e, respectively; Af denotes the number of excitable core electrons per unit volume. 
Only electric-dipole transitions will be retained between spin-orbit coupled inner orbitals, 
</5fj;-j^^. (r, s), localized around site i? = 0, and spin-polarized, spin-orbit coupled band states, 
tjjk (r, s) |13[. We consider only magnetic circular dichroism integrated over the two partners 




of a given spin-orbit split inner shell (nf) . For this, we have 



■AAcos0^(*o|afc4l*o) E E 



(1) 



Here, ep 



nl lii 



are the two threshold energies and Sc is a cut-off, positioned far above the 



top of the valence band so that nothing would change if £c were increased by the spin-orbit 
splitting, £f^ii^i — EfiH-i of the inner level. The superscripts ± identify circular polarisations, 
and QiM — eYij\/(f ) r is the electric dipole moment; 9 is the angle between photon wave vector 
and magnetization direction (orbital quantisation), which we take along the z axis; denotes 
the ground state of the system, and aj. is a fermionic creation operator for band states. 
To evaluate (ipk\QiAi\ VfOjfnj 



we use 



fnljr 



(2) 



for the inner orbitals, that is a two-component function with the same radial dependence, 
iPfii{r), for j = I ~ ^ and I + ^- Notice that this is an excellent approximation as the radial 
probability densities, 47rr^[/„(£, r)^ +g^{e, r)^], for k=1 and k= — l—l differ appreciably only 
close to the nucleus where they are small Q; here fK{e,r) and gi^{e,r) are the solutions of 
the radial Dirac equations with k—1, e=ef^ii_i and — 1 — 1, e=Sf^ii^i for the two partners, 
respectively. For ipk we use an expansion in spin functions times spherical harmonics centered 
at the absorption site 



V-fc (r, s) = ^ 



i,,kMr) Yim (f) Cm, (s) 



(3) 



The inner orbitals are so localized that, in the region relevant to the integral (^ipt IQiAf | fnTjin 
the self-consistent field for a band electron is dominated by the Hartree contribution, which 
is centrally symmetric and independent of fc. As a consequence, the coefRcients in expansion 



factorise into normalization constants, 



and radial functions, 4>i{r), which depend 



only on the magnitude of angular momentum about the absorption site. Regarding relativistic 
effects, the argument given above for neglecting the j — l± ^ dependence of the inner radial 
function holds also for the radial functions of the band states. The exchange-correlation po- 
tential does depend on the spin and possibly the orbital character but, in the inner region, it 
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amounts to m^- and possibly Zm-dependent shifts, which are small compared with the radial 
kinetic energy, Sk — v (r) — I {I + 1) /r^, in that region. Compared with that kinetic energy, 
the band-energy range, Sc — sp, is also small. 

We now write the difference between the squared matrix elements for right- and left- 
circularly polarized light as 

^ '^1M:W (V'fe IQimI y^firjrn, ) {vnVjni, \Q*Im\ V-fc ) , 

M=±l 

insert expressions (2) and (^), and apply the Wigner-Eckart theorem. Notice that, owing to 
the sum over j in (1), the resulting expression is spin independent. The angular part is then 
recoupled with use of the transformation p5|] 



2H^1 / ,^lm-l -M '^l'm':l-M' 



El -i^l" —m" +l — m —m" /~tl" m" } ^ ^ ^ 

l"m" ^ 

with I" — 0, 1, 2, corresponding to isotropic absorption, circular, and linear dichroism, respec- 
tively. Hence, we obtain the result 

/ [k+{u;) ~ K~icj)] dw = TT^e^AAcos^ ^ if^RiiT (*o |4| *o) , (4) 

with the radial integral over the inner region defined as 

Rm= Vni{r) r^i{r) r^dr^ , (5) 

and 'Pniif)'^ r^dr = 1. The operator is given by 

k va rris 

The individual normalisations of and (^) are irrelevant when simulating the integrated 
dichroism (jj). This is because normalization of the band states to unity in the solid merely 
fixes the normalization of the product u\^^ fc0i('') in (H)- The usefulness of atomic sum- 
rules, however, stems from a separation into an atomic factor which is independent of the 
magnetisation direction, and a remainder which, for circular dichroism, is approximately the 
ground-state orbital angular momentum of the excited atom. Suppose normalisations could be 
defined such that u\^^ were equal to the i? = component of the eigenvectors, u^imm fc' 
for the states 

RlmrUs 

in a representation of orthonormal orhitals. Then, from and (^), ^o) would be 

the expectation value of the orbital angular momentum in the ground state, and the atomic 
sum rule would also hold in the solid. 

Our use of LMTOs is motivated by the following features: they constitute a minimal basis 
whose orthonormal representation and spherical-harmonic expansions about neighboring sites 
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are well known; their use in practical computations is well established P-pl, even for systems 
with appreciable electronic correlations [p^, and for systems with spin-orbit coupling and 
spin-polarization [ pT| ; the simple formalism for the orthonormal set, which we shall use below, 
was recently re-derived without resort to the approximations of taking the interstitial kinetic 
energy equal to zero and of dividing space into atomic spheres P-pd|. Moreover, the LMTOs 
have recently been generalized to Nth-order MTOs spanning the states in a broad energy 
range, the accuracy and range increasing with N, for a fixed basis-set size | pl| . 

MTOs are derived from an MT-potential, i.e. a superposition of atom-centered spherically- 
symmetric potential wells with ranges limited to about 0.7 times the distance to the nearest 
neighbour jioj . In accordance with , we shall use MTOs derived from a spin- and orbital- 
independent potential: V (r) = ^j^ vr (|r — R|). To construct an LMTO, one first solves the 
appropriate radial Schrodinger or scalar-relativistic Dirac equation for various energies in the 
band region, thus obtaining the radial functions cfiRi (s^r). Each of these is then multiplied 
by Yim (?) and, if we are using the atomic-spheres approximation (ASA), truncated outside 
the atomic sphere. If not, they are augmented continuously with tails; these are localized 
(screened) solutions of the wave equation with the correct energies and are excluded from any 
mner region 1-0. As a result, we obtain the so-called truncated or kinked partial waves, 
0it/m(£,r - R). Now, the LMTO, 



XRlm (r - R) = (t)Rhn (r - R) + ^ <pR'l'm' (r - R') hR: I'm', Rim , (8) 

R'l'm 



centered at site R and with spherical-harmonic character Im, is defined as the corresponding 
partial wave, taken at energy Si, at the centre of interest, plus a smoothing cloud of the first 
energy derivatives, 

<i>R'i'm' (r-R') = d(j)R,i,m' (£,r-R') /ds\^^ , 

of partial waves at their own and at neighboring sites. (Here, and in the following, an omitted 
energy argument implies that e = e^.) In (|^), the expansion coefficients, h, form a Hermitian 
matrix which is approximately the band Hamiltonian with respect to e^, for the MT-potential 
used to generate the LMTO set. Specifically, since the LMTO is smooth, we may operate with 
(—A + V — Ev) on each term in (||) to obtain (—A + V ~ £)\4>{e)) = 0. Energy differentiation 
then yields (—A -\-V — e)\(l){e)) = |0(e)) and, as a result, 

(-A + V^(r) - £^)Xram(r - R) = ^ (I^R'l'm' -li!) hR'i'm' ,Rlm ■ (9) 

R'V m' 

As 9(/)(e, r)r(e) /de — (j){r)T + (r) T, with T—1, changing the energy-dependent normal- 
ization of a partial wave changes the shape of its energy-derivative function by adding some 
amount of 0(r) to it. This in turn changes the shape of the LMTOs via Eq. (H), but not the 
Hilbert space spanned by them. If each partial wave is normalized to one, we obtain a nearly 
orthonormal set since, in that case, the corresponding (pRimi^) is orthogonal to <j)Rim{'^)i as 
energy differentiation will reveal. Neglecting the overlap between partial waves at different 
sites (ASA), or using Lowdin orthogonalisation |p|-[lT||, one obtains 

{(t>Rlm I 4>R'l'm') — S RR.' 5mm' , {4>RlrrL \ 4'R'l'm''^ — 0. 

Insertion into ^ and (|^) finally shows that, in the nearly orthonormal representation, the 
LMTO Hamiltonian and the overlap matrices are given by {x\ — ^ + V — e^\x) = h and 
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(x I x) ~ 1 + hph^ respectively. (The off-diagonal elements of the matrix p \ 
neglected.) The truly orthonormal set is therefore 



<j)) h) (1 + hph) 2 



|X^) = + l-^^) = |X^) (1 + hph)-^ = + 
where the expansion matrix, 

= {1 + hph)~^ h {1 + hphyi , 

is the band Hamiltonian without spin polarization and spin-orbit coupling, and where 



may be 
(10) 

(11) 

-) = 



10) (1 + hph)^ and = ( " ^P) (1 + ^P^)^ ■ 

Next, we may work out the matrix elements of the exchange splitting and spin-orbit cou- 
pling in the orthonormal representation ([lO|), add them to e^^ -f /i^, and diagonalize to find the 
eigenvalues, Sfc, eigenvectors, j,, and band states (0). Expanding the latter in spheri- 

cal harmonics about the excited site using (|l^) , we are finally able to identify the coefficients 
■"Imm fe i'^ usual, an omitted subscript R imphes that R=0.) At first glance, it seems 

as if two radial integrals in (|^) are needed: one involving 0; (r), as contributed by the head of 
the LMTO, and the other involving (jji (r), as contributed mainly by the tails of neighboring 
LMTOs. However we observe that, when integrating the radial equation for the Z-channel 
outwards we may use the same initial condition for all energies. Hence, we obtain an energy- 
derivative function, (j)\{r), which is essentially excluded from the inner region, and whose 
contribution to the integral (||) may therefore be neglected . Since this procedure amounts 
to choosing a particular energy-dependent normalisation of the corresponding radial function, 
(j)\{s, r), the energy derivative function, (r) , must be a particular linear combination of the 
Yim projections of (j)im{v) and (pim (r) • These projections are independent of m in the ASA, 
where 4>im (^ii") = 0/ (e, r) (f), but only approximately independent when (e, r) is a 
Lowdin orthonormalized kinked partial wave. In the latter case, the m-dependence may be 
minimized through adjustment of the screening this dependence will be neglected in 

the present paper. Choosing to normalize (j)]^ (Sj r) to one at e^, we can express the linear 
combination which does not contribute to the radial integral (^) as a projection onto the 

orthonormal ^(/), (j)^ set 

-\<t>)o\ (12) 

Here, {(p \ </)') = o* is a matrix whose elements vanish unless R—R'=0 and I— I' = ^ ± 1, 
and whose off-diagonal elements are neglected together with any m-dependence. We thus 
eliminate 4>im (r) from ( p^ ) and find 



(1 + hphy 



(One should keep in mind that = \(j>) unless R=0 and I = l±l.) Identification of (|^) with 
(I) yields 

ul = {l- o'h){l + hphy^ui = (1 - o'h^ - ^h^ph^ + ...)ui. (13) 
and, hence, the final result for use in (ph is 





|2 _ 1 j_ 1 







2Re. 



. o; 



R'l'ra' 



Im.R'l'jn' ^Wl'rn'nisjk 
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0/ 



E 

R'l'm' 



(14) 



Re ■ 



hm, R'l'm' ^R'l'm'm^,k 



Rfl'm' 



The first term is contributed by LMTO heads only and gives the atomic smxi rule. Of the 
terms with R' ^ 0, the ones on the second line are LMTO head-tail contributions, and those 
on the remaining lines are tail-tail contributions. However, the sum of the terms on the 
second and third lines may also contribute to the atomic sum rule, as they depend on e^. 
[This dependence is cancelled by the dependence of the radial integral brought about by 
(t)i (£j/,r) in (||). Notice that to be able to neglect the m-dependence of the valence orbital 
in the radial integral, we have chosen this orbital as a partial wave rather than an LMTO, 
which has longer range.] To clarify this point, let us assume that the spin-orbit interaction is 
smaller than the exchange splitting and that the latter is fairly independent of k. In this case. 



E/ ^lm,R'l'm''^R'l'm'ms,k ^ {^k ^ v) ''^Imm ^ ,k ^ 



R'V m' 



where Cfc is the (doubly degenerate) band without spin-orbit and exchange couplings, and (| 
reduces to 



Imnis.k 



[l - {ck - e^)oJ]^ |wLm,,fc| 



I R'l'm' 



P)lmM'l'm' '^R'l'm',k 



We now realize that the deviation from the atomic sum rule is the contribution to the 
integrated dichroism (^ stemming from the second and further lines of (p^, after they have 
been minimized with respect to e^, that is, when is chosen as the centre of gravity of the 
unoccupied part of the /-projected density of band states 



Sul 



EkK 



, iV^^(v|/o|afc4|vI/o) J2 



This choice of Si, is the one which also minimizes the errors of the LMTO method. (When an 
/-independent is used, the o\Riii/{2l + l)-weighted average should be chosen). 

To summarise: to estimate the accuracy of atomic sum rules for x-ray dichroism in solids, 
we have examined the problem of x-ray absorption by band electrons, with emphasis on 
the interpretation of the total intensity of spectra obtainable with circular polarisation in 
magnetic systems. Using an orthonormal set of LMTOs, we have found corrections to the 
atomic results in the form of energy moments of the band. Applications of the approach to 
unpolarised and linear-dichroic spectra, together with a numerical determination of the actual 
size of the corrections in specific cases will be reported elsewhere. 



REFERENCES 



[1] A. F. Starace, Phys. Rev. B , 5 (1773) 1972; G. VAN der Laan and Thole B. T., Phys. Rev. 
Lett. , 60 (1977) 1988. 



R. Benoist, p. Carra and O. K. Andersen: Band structure and x-ray dichroism 7 



[2] B. T. Thole, Carra P., Sette F., and van der Laan G., Phys. Rev. Lett. , 68 (1943) 1992. 
[3] P. Carra, Thole B. T., Altarelli M., and Wang X. D., Phys. Rev. Lett., 70 (694) 1993. 
[4] P. Carra, Konig H., Thole B.T., Altarelli M., Physica, B 192 (182) 1993. 
[5] R. Q. Wu, Wang D., and Freeman A. J., Phys. Rev. Lett, 71 (3581) 1993; R. Q. Wu and 

Freeman A. J., Phys. Rev. Lett., 71 (1994) 1994. 
[6] G.Y. Quo, Ebert H., Temmerman W.M., and Durham P. J., Phys. Rev. B, 50 (3861) 1994. 
[7] See papers, for instance, by H. Ebert and by M.S.S. Brooks, and B. Johansson, in H. Ebert 

and SCHliTZ G. (Editor), Spin- Orhit- Influenced Spectroscopies of Magnetic Solids - Springer 

Lecture Notes in Physics, Vol. 466 (Springer Verlag, Berlin) 1996. 
[8] O. K. Andersen and Jepsen O., Phys. Rev. Lett., 53 (2571) 1994 

[9] O.K. Andersen, Jepsen O., and Krier G., in V. Kumar, Andersen O.K., and Mook- 

erjee a. (Editor), Lectures on Methods of Electronic Structure Calculations (World Scientific 

Publishing, Singapore) 1994, p. 63-124. 
[10] O. K. Andersen, Arcangeli C, Tank R. W., Saha-Dasgupta T., Krier G., Jepsen O., 

AND Dasgupta I., Mat. Res. Sac. Symp. Proc, 491 (3) 1998 
[11] O.K. Andersen, Saha-Dasgupta T., Tank R. W., Arcangeli C, Jepsen O., and Krier 

G., in H. Dreysse (Editor), Ab Initio Calculation of the Physical Properties of Solids - Springer 

Lecture Notes in Physics (World Scientific Publishing, Singapore) 535 
[12] The absorption coefficient is defined by ^^{uj) — {2uj / c) (uj) . 

[13] Translational symmetry is no prerequisite for our formalism. For states with translational sym- 
metry, fc is a joint label for the band index and the Bloch vector. 

[14] This approximation causes errors less than 1% for the L2,3 edges of 3d transition metals and 
the M4,5 edges of rare earths 

[15] D. a. Varshalovich, Moskalev A. N., and Kersonskii V. K., Quantum Theory of Angular 
Momentum (World Scientific Publishing, Singapore) 1988, p. 453. 

[16] V.I. Anisimov, Zaanen J., and Andersen O.K., Phys. Rev. B, 943 (3) 1991; A.I. Liecht- 
enstein, Anisimov V.I., AND Zaanen J., Phys. Rev. B, 52 (R5467) 1995. 

[17] H. Ebert, Phys. Rev. B, 38 (9390) 1988; I.V. Solovyev, Liechtenstein A.I., Gubanov 
V.A., Antropov V.P., AND Andersen O.K., Phys. Rev. B, 14 (414) 1991. 

[18] O.K. Andersen, Phys. Rev. B, 12 (3060) 1975; 



